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We perform QCD fit of the nucleon strange and anti-strange sea distributions to the 
neutrino and anti-neutrino dimuon data by the CCFR and NuTeV collaborations, sup- 
plemented by the inclusive charged lepton-nucleon Deep Inelastic Scattering (DIS) and 
Drell-Yan data. The effective semi-leptonic charmed-hadron branching ratio is con- 
strained from the inclusive charmed hadron measurements performed by the FNAL- 
E531 and CHORUS neutrino emulsion experiments as = (8.8 ± 0.5)%. We obtain 
a strange sea suppression factor k(20 GeV^) = 0.62 ± 0.04(exp.) ± 0.03(QCD). An x- 
distribution of total strange sea obtained in the fit is slightly softer than the non-strange 
sea, and an asymmetry between strange and anti-strange quark distributions is consis- 
tent with zero (integrated over x it is equal to 0.0013 ± 0.0009(exp.) ± 0.0002(QCD) at 
the scale of 20 GeV^). 

At Bjorkcn variable x < 0.2 the strange quarks give important contribution to the total 
quark sea in the nucleon. Therefore accurate determination of the strange sea is necessary 
for interpretation of the precise hadron-coUider and fixed-target experimental data. In par- 
ticular, the presence of a small positive s — s asymmetry in the nucleon may explain [2] 
recent anomalous result on the weak mixing angle in the NuTeV experiment [3]. The 
best constraint on the strange sea comes from the DIS neutrino-nucleon dimuon produc- 
tion. This process stems from the charged-current production of the charm quark decaying 
semileptonically with a secondary muon in the final state. This mechanism is particularly 
sensitive to the strange quark distributions, as the contribution from non-strange quarks 
are greatly suppressed due to smallness of the corresponding Cabibbo-Kobayashi-Maskawa 
(CKM) matrix elements. This allows to disentangle the strange distribution from other 
quarks in the global PDF fit, which includes the dimuon data. 

We extract the strange distribution from the dimuon data by the CCFR and NuTeV 
experiments [4j [5l [6] . Those data were collected in the runs with the neutrino and anti- 
neutrino beams and allow separate determination of the strange and anti-strange distribu- 
tions at 0.01 < X < 0.3 and the values of the momentum transfer 1 GeV^ ^ ^ 200 GeV^. 
In order to consistently account for the contribution from other quarks and giuons, the 
dimuon data are supplemented by the inclusive DIS data and the fixed-target Drell-Yan 
data (cf. Ref . [7] for details of the data selection) . The analysis is performed with account 
of the NNLO correction to the parton evolution and the massless coefficient functions [9] . 
The heavy-quark DIS production is calculated in the 3-flavour factorization scheme with the 
perturbative corrections up to the NLO [10 . The latter are illustrated in Figure[T]for the rep- 
resentative kinematics of the data used in the fit and the factorization scale ji ~ + 
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y=0.5, Q^=10 GeV^ 




Figure 1: The values of the corrections on the calculated dimuon cross sections employed in our 
analysis (solid lines; the NLO QCD corrections, dotted-dashes: the nuclear correction for the case of 
iron target) for the neutrino (left panel) and anti-neutrino (right panel) beam at the representative 
kinematics of the NuTeV and CCFR dimuon data, = 10 GeV^ y = 0.5. 



where mc is the charm quark mass. At small x the NLO corrections of Ref. |10| reduce the 
value of calculated cross sections therefore they increase the strange distributions extracted 
from the data. Corrections for nuclear effects in the neutrino-nucleon DfS, which include 
the Fermi motion and binding, neutron excess and shadowing, nuclear pion excess and the 
off-shell effects [TT] are also taken into account in our analysis. The value of nuclear correc- 
tions depends on the beam type (cf. Figure [T]) , therefore they affect the charge asymmetry 
in the strange distribution x\s(x) — Other details of the analysis including treatment 

of the data and their uncertainties can be found elsewhere [T^] . 

Assuming universality of the charmed hadron fragmentation function, the dimuon pro- 
duction cross section is related to the charm production cross section as follows 

dxdy dxdy ' 

where y is inelastisity, ry^ is the acceptance correction accounting for the cut iJ^ > E^ im- 
posed on the muon energy to reject background, and is the effective semileptonic charm 
quark branching ratio. The value of = ^ f\iBr{h — » ^X), where /h and Br{h — > /^^) 
are the rate and the semileptonic branching ratio of the charmed hadron h produced in 
the (anti)neutrino-nucleon scattering, respectively, h = , , D'^ , . The value of 
depends on the beam type and energy, due to f^. For kinematics of the NuTeV and 
CCFR experiments a value of B^ = (9.19 ± 0.94)% averaged over neutrino and antineu- 
trino beams was obtained in Ref. |13j using the charmed hadron fractions measured by the 
Fermilab-E-531 experiment. In our analysis we fit the value of i?^ simultaneously with 
the PDF parameters. In this case the strange quark distributions are constrained by the 
Q^-slope of the dimuon cross sections, which depend on the strange sea magnitude due 
to the PDF evolution, and the value of B^ is constrained by the magnitude of the cross 
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Figure 2: The ilcr uncertainties in the C-even (left panel) and C-odd (right panel) combination of 
the strange sea distributions determined in our fit with (solid lines) and without (dashed lines) the 
constraint on at the value of factorization scale /i = 3 GeV. The difference between distributions 
obtained in these two variants of the fit is shown by the dashed-dotted lines. 



sections. This approach provides a consistent determination of and allows to check 
the beam- and energy-dependence of i?^ for the kinematics of experiments used in the fit. 

We do not observe statis- 
tically significant beam- 
and energy-dependence of 
Bp. The averaged over 
beam type and energy 
value of = (9.1 ± 
1.0) % obtained in our fit 
is in good agreement with 
the result of Ref. [13]. 
The value of is anti- 
correlated with the mag- 
nitude of the strange sea, 
therefore an additional constraint on i?^ can improve the strange sea determination. The 
value of Bf^ obtained in Ref. [T^ does not provide efficient constraint since its uncertainty is 
consistent with the corresponding uncertainty in our fit. The uncertainty in i?^ can be re- 
duced by utilizing the results of the recent measurement of the charmed hadron fractions by 
the CHORUS experiment [111 [15] as well as the updated determination of the semileptonic 
branching ratios Ref. [16]. The value of i?^ obtained with these inputs is given in Table 1. 
The uncertainty of this determination is about factor of 2 smaller than that of Ref. [13] and 
the central value of B^ is somewhat lower (cf. Table 1.). In the variant of our fit with the 
constraint of Table 1 for the beam energy cut > 30 GeV, the central values of the strange 
and antistrange distributions are smaller, correspondingly. In this fit the uncertainty in the 
C-even combination x[s{x) + s{x)]/2 is about factor of 2 smaller than for the fit with no 
constraint on B^ (cf. Figure[2]). For the C-odd combination x[s{x) — s{x)], the effect of the 



Measurement 


E,, > 5 GeV 


> 30 GeV 


E531 charmed fractions [ISJ 


7.86 ± 0.49 


8.86 ±0.57 


CHORUS direct [H] 


7.30 ± 0.82 


8.50 ± 1.08 


CHORUS charmed fractions [15] 


9.11 ±0.93 




Weighted average 


7.94 ± 0.38 


8.78 ±0.50 



Table 1: Semileptonic branching ratio _B^(%) from direct measure- 
ments in the E531 and CHORUS emulsion experiments. The last row 
corresponds to our weighted average. 
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constrained is weaker because of a partial cancellation in the difference. 

In the fit with the constrained B^, we obtain the values of /NDP of 63/89 and 38/89 
for the CCFR and NuTcV data sets, respectively. This is statistically consistent with the 
effective number of degrees of freedom for these experiments introduced in Ref. [51 , in order 
to take into account statistical correlation between the different data points that is roughly 
twice smaller than NDP. The value of the strange sea suppression factor k, which is the 
ratio of sum of momentum carried by the strange and anti-strange quarks to the one carried 
by the sea up- and down-quarks, is given in Figure [3l Due to the QCD evolution the value 
of K, rises with the factorization scale /i. The value of k(20 GeV^) = 0.48 ^[J'^g was obtained 
in the NLO fit of Rcf. ^ based on the CCFR dimuon data. In our fit we get a bigger value, 
k(20 GcV^) — 0.62 ± 0.04. This is because the non-strange sea distribution used in Rcf. [3| 
is not consistent with the distribution of the present fit as well as with other modern PDF 
sets |171 [T8] . Our value of k, is somewhat smaller than the value for the MSTW08 PDFs 
set of Ref. [17] and is somewhat bigger than one for the CTEQ6 PDFs set of Ref. [18]. 
In both cases, however, the discrepancy is within uncertainty in our determination of k, 
which is obtained from propagation of the statistical and systematical errors in the data 
with the account of correlations. As can be seen from Figure [3l the C-even strange sea is 
somewhat steeper than the non-strange sea and at small x overshoots the latter. At a; < 0.2 
the strange sea is enhanced due to the NLO corrections to the charm production coefficient 
functions. Were these corrections dropped, we obtain the value of k(20 GeV^) = 0.55±0.13. 
The uncertainty in k(20 GcV^) due to the higher-order QCD corrections, estimated as the 
change due to variation of the dimuon cross section factorization scale from y'Q^ + mf to 
Q, is 0.03. This is comparable to the experimental uncertainty in k. In this sense the 
theoretical and experimental uncertainties in our determination of k are balanced, and, 
in order to further improve it, more experimental data arc needed and the NNLO QCD 
corrections have to be taken into account. 

In a variant of our fit with only the NuTcV dimuon data, the strangeness charge asymme- 
try is somewhat positive. This is consistent with the results of the NLO analysis of Rcf. [6]. 
The CCFR data prefer somewhat negative asymmetry and the result averaged over these 
two experiments is consistent with in a wide range of ^, even despite it rises with /i due 
to the NNLO corrections to the PDF evolution [l9j. Furthermore, if we fix the strange sea 
asymmetry at 0, the value of rises by about 1 unit only. For the fit based on the com- 
bined NuTeV and CCFR dimuon data we obtain S- = x[s{x) -l{x)]dx = 0.0013±0.0009 
at /z^ = 20 GeV^, that is also consistent with within the uncertainties. The theoretical 
uncertainty in due to the factorization scale variation is 0.0002, much smaller than both 
the experimental uncertainty and the theoretical uncertainty in k. 

In summary, we obtain the strange sea suppression factor k(20 GeV^) = 0.62±0.04(exp.)± 
0.03(QCD), the most precise value available. The ^-distribution of the total strange sea is 
slightly softer than the non-strange one, and the integral strange sea charge asymmetry 
S'-(20 GcV^) = 0.0013 ±0.0009(exp.)±0.0002(QCD) is consistent with within uncertain- 
ties. 
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Figure 3: Left panel: The ±lcr band for the strange sea suppression factor k with respect to 
the factorization scale jj, as obtained in our analysis (solid lines) in comparison with ones for the 
MSTW08 (dashes) and CTEQ6 (dots) PDF sets. Right panel: The ±lcr band for the C-even 
combination of the strange sea distributions determined in our fit (solid lines) compared to the 
non-strange one scaled by k (dashes) at = 3 GeV. 
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